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ABSTRACT 

We present estimates of CN and Mg overabundances with respect to Fe for 
early-type galaxies in 8 clusters over a range of richness and morphology. Spectra 
were taken from the Sloan Digital Sky Survey (SDSS) DR1, and from WHT 
and CAHA observations. Abundances were derived from absorption lines and 
single burst population models, by comparing galaxy spectra with appropriately 
broadened synthetic model spectra. We detect correlations between [Mg/CN] and 
[CN/Fe] and cluster X-ray luminosity. No correlation is observed for [Mg/Fe]. 
We also see a clear trend with the richness and morphology of the clusters. This is 
interpreted given varying formation timescales for CN, Mg and Fe, and a varying 
star formation history in early- type galaxies as a function of their environment: 
intermediate-mass early-type galaxies in more massive clusters are assembled on 
shorter timescales than in less massive clusters, with an upper limit of ~ 1 Gyr. 

Subject headings: cosmology: observations — galaxies: abundances — galaxies: 
clusters: general — galaxies: formation — galaxies: stellar content — X-rays: 
galaxies: clusters 

1. Introduction 

A key question for scenarios of galaxy formation is whether galaxies formed in single 
structural "monolithic" events (Larson 1974) or by a series of "hierarchical" processes (Press 
& Schechter 1974; White & Frenk 1991) in which large galaxies built up from smaller ones. 
Many structural and dynamical properties of galaxies in clusters are explained in this sce- 
nario, though problems remain: the absence of the predicted mass cusps in the centres of 
ellipticals and bulges, and the prediction of far more satellite galaxies than those observed. 
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Stellar populations offer a fossil record of the formation and evolution of galaxies, most 
clearly in elliptical galaxies, and stellar population studies provide very strong constraints 
on the principal galaxy formation scenarios. It is hard to reconcile the hierarchical models 
with the result that massive galaxies show significantly larger mean luminosity weighted ages 
than their smaller counterparts (Kauffmann et al. 2003). 

Understanding stellar populations in early-type galaxies as a function of the environment 
can provide answers to the puzzle. The present observational base is small. Only three 
clusters have been observed for a detailed stellar populations analysis: Virgo (e.g. Vazdekis et 
al. 2001a), Coma (e.g. J0rgensen 1999) and Fornax (Kuntschner & Davies 1998). Extending 
the study to a large number of clusters covering a range of richness and morphology is 
mandatory. 

Past studies of clusters used the original Lick/IDS spectral indices (Worthey 1994) whose 
reliability is limited by their resolution dependence, as uncertain corrections for broadening 
and instrumental effects are needed. It is better to use modelled integrated spectra, which 
can be broadened to match galaxy velocity dispersion, a. Such models have been developed 
by Vazdekis (1999) and allow an accurate separation of age and metallicity, and subsequent 
individual abundance derivations. 

The study of the element abundance ratios in elliptical galaxies within distinct clusters 
should be a powerful discriminant between different star formation histories (e.g. Worthey 
1998). In particular, overabundances of [Mg/Fe] compared with the solar ratio have been 
found in massive elliptical galaxies (Peletier 1989; Worthey et al. 1992; Vazdekis et al. 1997). 
These have been interpreted via several possible scenarios based on the fact that Mg is mainly 
produced in Type II supernovae (Faber, Worthey, & Gonzalez 1992; Matteucci 1994), and 
include different star formation rates (SFR) and a time dependent IMF. 

Differences in the abundances of C and N as a function of the environment have been 
recently suggested by Sanchez-Blazquez et al. (2003), who found striking spectral differences 
between field elliptical galaxies and their counterparts in the central region of the Coma clus- 
ter. Galaxies in the denser environment showed significantly lower CN 2 and C4668 absorption 
strength. Here we explore these differences by extending the study to a larger number of 
clusters over a range of richness and morphology, applying the new analysis techniques to 
derive abundance ratios. We have assumed a flat Universe with H = 75 km s _1 Mpc -1 and 
g = 0.5. 
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2. Data 

Sloan Digital Sky Survey (Stoughton et al. 2002) spectra were obtained using a multi- 
object, 3" diameter fiber spectrograph. Exposures ranged from 45 to 75 minutes. All the 
data processing was performed with automated SDSS software. Redshifts were measured 
on the reduced spectra by an automated system, which models each galaxy spectrum as a 
linear combination of stellar populations. We measured independently the redshifts of the 
galaxies used in this study by crosscorrelating each galaxy spectrum with our SSP synthetic 
model spectra. We found no significant differences between SDSS redshift values and ours. 

From the SDSS Data Release 1 database we selected galaxy spectra according to the 
following criteria. They must: 

i) Belong to an Abell cluster. We included this criterion because the richness and the 
morphology of Abell clusters are uniformly defined and described in the literature. Also, 
X-ray luminosity values are available. 

ii) Belong to the early-type galaxies catalogue of Bernardi et al. (2003). This sample 
has ~ 9000 early-type galaxies, in the redshift range 0.01 < z < 0.3, selected from the 
SDSS spectroscopic database using morphological and spectral criteria. The mean spectrum 
signal-to-noise per pixel is 16. 

Hi) Have S/N per pixel greater than 15. 

iv) Have velocity dispersion in the range 150 km < a < 250 km s" 1 . Galaxies outside 
this range of a were rejected because of the completeness of the sample: not all clusters had 
spectra of dwarf and/or giant elliptical galaxies because of the inner limitation of SDSS data 
(for dwarfs) and the morphology of the clusters (for giants). Also, the quality of the spectra 
of the faintest galaxies was too low for our analysis requirements, a values were obtained 
from Bernardi et al. (2003). 

Using these criteria, we obtained a total of 55 galaxies distributed in 6 clusters. The 
clusters are: A257, A279, A655, A1238, A1650 and A2050. Their redshift values vary in the 
range 0.07 < z < 0.13 and they cover a range of richness and morphology. See details in 
Table 1. 

For comparison and completeness, we added to our SDSS data high-quality long-slit 
spectra of early-type galaxies in Coma and Virgo clusters (for details, see Sanchez-Blazquez 
et al. 2003). To compare them with those of the SDSS sample, we extracted spectra along 
the slit simulating a circular aperture (distance weighted coadded spectra) of radius 1.5" at 
z = 0.1. This aperture translates into apertures of radius 6" for Coma and 37" for Virgo. 
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3. Galaxy measurements and results 

To derive mean luminosity-weighted ages and metallicities, we compared selected ab- 
sorption line strengths with those predicted by the model of Vazdekis (1999). This model 
provides flux-calibrated spectra in the optical range at a resolution of 1.8 A (FWHM) for 
single-burst stellar populations. This way, we can transform synthetic spectra to the resolu- 
tion and dispersion of the galaxy spectra instead of the opposite, as required while working 
in the Lick system. Selected absorption indices were CN 2 , Mg 2 (Worthey et al. 1994) and 
Fe2 (defined as Fe2 = Fc4383 + Fc527 ° ) _ \y e used these features because of their low sensitiv- 
ity to variations in S/N (Cardiel et al. 2003) and velocity dispersion (we have estimated 
A (index) / index < 0.15, for A(er) = 300 km s -1 ). This way we avoid possible variations in 
the index value, as a may vary as a function of r, due to the fact that SDSS spectra provide 
the light integrated within the fibers of the spectrograph. 

Plots of the strengths of the selected indices versus H/5 provide close to orthogonal model 
grids, allowing us to estimate accurately galaxy mean ages as well as relative abundances of 
the different elements. Figure 1 illustrates this method for the galaxies of clusters A1238 and 
A655 (two clusters with extreme values of X-ray luminosity). We will refer to the metallicities 
derived in the diagrams CN 2 -H/3, Mg 2 -H/3 and Fe2-H/3 as Z C n, Z Mg and Z Fe , respectively. 
Since the CN 2 index is strongly dominated by C and N, the Mg 2 index is governed by Mg 
and the Fe2 index by Fe (Tripicco & Bell 1995), these metallicities must be close to the 
[CN/H], [Mg/H] and [Fe/H] abundances, and [Zc^/Z Fe \, [Z Mg /Z Fe ] and [Z Mg /Z CN ] are then 
estimates of the abundance ratios [CN/Fe], [Mg/Fe] and [Mg/CN] for each galaxy. Note that 
an extrapolation of the model grids is required for some galaxies to obtain the abundances 
of CN and Mg, since the models extend only to [M/H] = 0.2. It is worth noting that certain 
galaxies fall below the model grids, which can be attributed to the fact that absolute age 
determination is subject to model uncertainties (see Vazdekis et al. 2001b; Schiavon et al. 
2002). We neglect the possible effect of nebular emission on Hg since no [O III] A5007 emission 
is detected in our galaxy spectra. Nevertheless, assuming an upper EL3 emission correction of 
~ 0.5 A (Kuntschner et al. 2002), although this would give rise to a significant reduction in 
the mean age of the stellar populations of the oldest galaxies, the net effect on the abundance 
ratios would be no more than ~ 0.05 dex, for the most affected galaxies in a cluster. Note 
that this correction for the B.(3 index is larger than the one obtained by varying the model 
prescriptions, e.g. with a-enhanced isochrones and atomic diffusion included (Vazdekis et 
al. 2001b), for inferring ages in better agreement with the current age of the Universe, for 
the oldest galaxies in our sample. However we do recognize that the ages quoted here, and 
in current articles dealing with stellar population modelling, may well span an older range 
than the true age range of the populations observed. This would indicate some flaw or flaws 
in the present quantitaive understanding of stellar evolution, and is a problem recognized by 
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the community working in population synthesis. 

Figure 2 shows the measured abundance ratios for the galaxies in the whole sample of 
clusters as a function of the velocity dispersion. To derive representative relative abundances 
ratios for each cluster, and since relative abundances correlate with a, we have fitted to each 
cluster a straight line with a fixed slope and a varying intercept. We have assumed that this 
slope corresponds to that found for the Coma cluster in the considered range of a. Note 
that we do not have any, statistically significant, evidence of a variation of the slope between 
the different clusters. In Table 1 we list the derived relative abundances at a fixed velocity 
dispersion of a = 200 km s" 1 for each cluster. 

As a quantitative indicator of the masses of the clusters, we have used X-ray luminosities. 
Lx values were taken from Ebeling et al. (1998) and Ledlow et al. (2003), adapted to the 
cosmological parameters used by the latter. 

Figure 3 shows the values of [Zcn/^Fc], [^Mg/^Fe] and [Zy[g/Zc^} versus X-ray luminos- 
ity for each cluster. We find clear correlations between [CN/Fe] and [Mg/CN] values and 
X-ray luminosity, with probabilities of no correlation of 0.046 and 0.002, respectively. On 
the contrary, no correlation is found for [Mg/Fe]. It is worth noting that the significant point 
in these relations is the relative differences in abundance ratios, not their absolute values. 
It is noticeable that the abundance ratio values also correlate with the richness class and 
the morphological type (see Table 1), with probabilities of no correlation of 0.05 and 0.02, 
respectively. 

4. Discussion 

The correlations can be interpreted in terms of the different formation timescales for 
each element, and the different star formation histories of early-type galaxies as a function 
of their environment. 

Magnesium is ejected into the interstellar medium (ISM) by Type II supernovae (SNe II) 
on short timescales (< 10 Myr). On the other hand, the iron-peak elements are the products 
of SNe la, which occur on timescales of ~ 1 Gyr. Between the two extremes, although there 
are recent suggestions that most of the C come from massive stars (Akerman et al. 2004), C 
and N are mainly ejected into the ISM by low- and intermediate-mass stars (Renzini & Voli 
1981; Chiappini, Romano & Matteucci 2003), leading to CN formation on timescales longer 
than for Mg but shorter than for Fe. Furthermore, several authors (e.g. Ellis et al. 1997; 
Stanford, Eisenhardt, & Dickinson 1998) argue that early-type galaxies are old and passively 
evolving systems. In any case, the luminosity-weighted ages derived from our model grids 
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confirm that the galaxies are significantly older than the formation timescales of the different 
species. So, if we find substantial differences in the abundance ratios of these elements which 
depend on the physical properties of the environment, these must be due to the fact that 
galaxies are assembled on different timescales as a function of their environment. 

In this framework, the constancy of the [Mg/Fe] values is explained in terms of the great 
difference in the formation timescales of the two elements: the galaxies are fully assembled 
before Type la SNe can significantly pollute with Fe the ISM of the smaller galaxies before 
merging, and right after Mg is fully ejected. Since [Mg/Fe] is found to be constant with the 
X-ray luminosity of the clusters (see Fig. 3), which is an indicator of their mass, we conclude 
that this ratio is independent of the environment. Similar results for the [Mg/Fe] ratio have 
been obtained by other authors (J0rgensen 1999; Kuntschner et al. 2002; Sanchez-Blazquez 
et al. 2003), by studying field and Coma cluster elliptical galaxies. 

However, when considering species with less disparate formation timescales, such as 
CN and Fe, or CN and Mg, clear correlations are found between abundance ratios and the 
environment, as shown in Fig. 3. The fact that [CN/Fe] decreases with the cluster X-ray 
luminosity, and that [Mg/CN] increases with it, suggests that galaxies in more massive 
clusters are fully assembled on shorter timescales than those in less massive clusters. We 
show that this difference is large enough to produce measurable variations of the abundance 
ratios of galaxies in more or less massive clusters. 

The result that there exist relative differences in the assembly timescales of the galaxies 
due to the properties of the environment is qualitatively in agreement with the hierarchi- 
cal models. Discrepancies appear, however, when considering the absolute values of such 
timescales. The fact that [CN/Fe] abundance ratio is not constant implies that early-type 
galaxies are fully assembled on timescales around the massive release of CN into the ISM. 
Hierarchical models, on the contrary, predict longer assembly timescales. 

Other scenarios have been explored in order to explain the differences in abundance ratio 
values as functions of the environment (see Sanchez-Blazquez et al. 2003). These include a 
decrease in the stellar giant/dwarf ratio in high-density environments, with respect to low- 
density ones, which would lead to lower index values in the latter. But model calculations 
have shown that the differences due to a variation in the IMF are too small to produce 
the observed variations. Also, a difference in the luminosity-weighted mean age between 
high- and low-density environments has been proposed. Models show that, to account for 
the differences, the galaxies in high-density environments must be ~ 8 Gyr younger than in 
low-density ones. This contradicts previous studies that suggest that galaxies in high-density 
environments are, in any case, older than those in low-density environments (Kuntschner et 
al. 2002). So our interpretation here seems the most consistent. 
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It is noteworthy that we have found a dependence of the abundance ratios with several 
properties of the environment, both quantitative (i.e. X-ray luminosity values) and qualita- 
tive (i.e. richness and morphological types) as shown in Table 1. This lends strong support 
to the basic hypothesis that the characteristics of the environment affect the evolution of the 
galaxies. The relations we have found set clear constraints on models of chemical evolution 
and galaxy formation. 

The authors thank A. Aguerri, X. Barcons, L. Carigi, C. Gutierrez, R. Peletier and 
the anonymous referee for useful comments. We acknowledge grant AYA2001-0435 from the 
Spanish Ministry of Science and Technology. The SDSS is managed by the ARC for the 
Participating Institutions. 
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Fig. 1. — Model grids for the galaxies of clusters A1238 (filled circles) and A655 (open 
circles). These two clusters have extreme values of X-ray luminosity. From left to right: 
Plots of H/5 vs. the metallicity indices CN 2 , Mg 2 and Fe2 (defined in the text). The galaxy 
velocity dispersions increase from top to bottom, in the three bins of a quoted in the boxes, 
to which galaxy and model spectra have been broadened. Overplotted are the models by 
Vazdekis (1999). Lines of constant [Fe/H] = —0.7, —0.4, 0.0 and 0.2 are shown by dot-dashed, 
dotted, solid and dashed lines, respectively. Thin dotted horizontal lines represent models of 
constant age, quoted in gigayears. The ages become increasingly uncertain for values greater 
than 8 Gyr (see Worthey 1994; Vazdekis 1999) but these uncertainties have no significant 
effects on relative abundance determination. 
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Fig. 2. — Abundance ratios [Zcn/^Fc] (left), [Zyig/Z-pe] (centre) and [ZMg/ZcN] (right) versus 
galaxy velocity dispersion, a. Each point corresponds to one galaxy within a cluster. Errors 
are computed from the measured index error bars shown in Fig. 1. Galaxy clusters are 
ordered by decreasing X-ray luminosity, from top to bottom. The dotted straight lines show 
the level of relative abundances for each cluster. See the text for more details. 
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Fig. 3. — Cluster X-ray luminosity versus overabundance values of [ZcN/-^Fe] (top), [ZMg/Zp e } 
(middle) and [Zu g /Zc^] (bottom). Circles, squares and triangles indicate clusters with rich- 
ness classes of 1, 2 and 3, respectively. Each point corresponds to one individual cluster and 
represents the interpolated abundance ratio for a = 200 km s^ 1 . The Spearman rank-order 
correlation coefficients and its significance values are written in top and bottom panels. 
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Table 1. Properties and measurements of analyzed clusters 



Cluster 


z 


Rich. a 


Morph. b 




[Zcn/ Zpe] 


[Z Mg /Z Fe ] 


A0279 


0.080 


1 


II III 


0.08 


0.43 


0.33 


A1238 


0.073 


1 


III 


0.15 


0.61 


0.44 


Virgo 


0.004 


1 


III 


0.30 


0.48 


0.36 


A0257 


0.070 


1 


II III 


0.31 


0.45 


0.35 


A2050 


0.118 


1 


II III 


1.22 


0.40 


0.31 


Coma 


0.023 


2 


II 


1.80 


0.39 


0.38 


A0655 


0.127 


3 


I II 


1.97 


0.30 


0.43 


A1650 


0.085 


2 


I II 


2.01 


0.40 


0.39 



a Abell richness class. b Bautz-Morgan morphological type. c X-ray luminosity in units of 
10 44 erg s" 1 . 



